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Abstract 

With the completion of a 24 operational satellite constellation , GPS is fast approaching the 
critical milestone f Full Operational Capability (FOC). Although GPS is well capable of providing 
the timing accuracy and stability figures required by system specifications , the GPS community will 
continue to strive for further improvements in performance. 

The GPS Master Control Station (MCS) recently demonstrated that timing improvements are 
always composite Clock , and hence , Kalman Filter state estimation f providing a small improvement 
to user accuracy. 


INTRODUCTION 

Though well capable of meeting and/or exceeding customer expectations, the GPS Master 
Control Station (MCS) will continuously search for safe and efficient methods for improving 
GPS timing accuracy and stability performance. The most recent improvements have focused 
on fine tuning the Continuous Time Update Process Noises (a.k.a. q s) for all GPS satellite 
frequency standards. 

Process noises are nothing new to the timing community. Many time scale algorithms update 
these parameters dynamically for their respective systems. As in many Kalman Filters, the 
Defense Mapping Agency (DMA) periodically reviews their q values for their OMNIS compu- 
tation program. OMNIS, like the MCS Kalman filter, estimates the ephemeris, solar, and clock 
states for 25 GPS satellites!^. However, up until 6 Oct 94, the timing community had never 
undertaken the task of re-qing an entire operational GPS constellation in the MCS Kalman 
Filter. 

Thanks to the generous input from several outside agencies, we now employ process noise 
values that are unique to the individual characteristics of the 25 operating frequency standards 
on orbit. Perhaps more importantly, we now also have the precise data, know-how, tools, and 
procedures to safely and efficiently review and update our q values on a periodic basis. 
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RUBIDIUM CLOCK ESTIMATION 


Each GPS satellite uses one of two different types of atomic clocks to provide a stable output 
frequency, to, in turn, generate accurate navigation signals. The majority of Block II/IIA GPS 
satellites currently use one of two available Cesium (Cs) frequency standards. Orbiting Cs 
clocks demonstrate reliable performance, with one— day stabilities ranging between 0.8 x 10 13 
to 2.0 x 10 — 13 I 13 * 14 » 1S I. The drift rate term for a Cs frequency standard is typically on the order 
of 1 x 10 -20 s/s 2 or less. Such a small drift rate term, an order of magnitude smaller than our 
time steering magnitude, has negligible effects on GPS timing (hundredths of a nanosecond 
over one day). Because of its relatively insignificant effect on frequency estimation, the MCS 
currently fixes the drift rate estimate to zero for all Cs frequency standards (on-orbit and 
ground based). 

Two Rubidium (Rb) clocks also reside on each Block II/IIA satellite. Rb clocks do exhibit a 
significant aging characteristic, typically on the order of 1 x ICC 18 s/s 2 . However, if a filter 
properly corrects for drift rate, the typical one-day frequency stability of a Rb clock state is 
significantly better than that of a Cs (0.6 x 10~ 13 versus 1.0 x 10 -13 ) H3,14,151 Unfortunately, in 
the past, our Kalman Filter had difficulty estimating drift rate. As a result, Rb clock estimates 
have had somewhat large variances, causing, in turn, increased difficulty in estimating frequency. 
Although a Rb clock itself is usually more stable than a Cs at one day, the stabilities of the 
MCS’s Kalman filter Rb clock states have, in the past, been worse than those for Cs clocks. 

This filter instability has impeded the MCS from incorporating their inherently better stability 
into GPS time calculations. Consequently, the timing community has been uneasy about using 
Rb clocks in GPS. Of the first 24 operational satellites, we initialized only three with Rubidium 
clocks. 

Despite this reluctant attitude towards using Rubidium clocks, many have realized that as 
Cesium clocks reach their respective ends of operational life, we will have no choice but to 
use more Rubidium clocks. In any case, it seemed counterintuitive that GPS was not making 
the most use of our most stable clocks. In early 1994, the 2 SOPS Navigation Analysis Section 
began tackling this long-standing concern. Because the problem resided in estimation, as 
opposed to physical clock performance, the Kalman Filter really only needed a fine tuning. 


Deriving New Rubidium Clock qs 

The MCS Kalman filter performs recursive time and measurement updates of the state residuals 
and covariances. In pure prediction, the clock state covariances are functions of the system qs 

H*l : 


P = 


q\t + q 2 t 2 + q^t? / 20 
<72* 2 /2 + </3* 4 /^ 

q 3* 3 /6 


q 2 t 2 /2 + q 3 t 4 /8 

<72* + </3t 3 /3 
<73* 2 /2 


<73* 3 /6 
<73* 2 /2 
<73* 


( 1 ) 


The Naval Research Laboratory (NRL) produced a report for 2 SOPS (ALL-5, 27 Jan 94), on 
SVN25. The report included a series of drift rate plots for the Rb clock that was active from 
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Mar 1992 until Dec 1993. NRL plotted 5, 10, 20, and 30-day averaged values for drift rate! 11 !. 
In analyzing the 30-day average plot, we noticed that the drift rate changed significantly more 
during the first 90 days than during the remaining operational time [Figure 1], 

From the above P matrix, in pure filter prediction, the system variance for drift rate is the 
scalar time product of q 3 : 

C 3 = t l3( r )(^) ( 2 ) 

Using the above equation, along with the NRL data, we derived new q values, both from the 
90-day initialization period, and from the remaining period, and we compared these to the old 
system q values: 


q value OLD INITIAL NEW NORMAL 

Drift Rate (</ 3 ) 9.00 x 10- 42 s 2 /s 4 1.35 x 10- 43 s 2 /s 5 6.66 x 10 -45 s 2 /s 5 

We also looked at calculating a new drift (frequency) q value. The old Rb q value for drift, 
4.44 x 10 _32 s 2 /s 3 , was the same as that for Cs. We chose </2 = 2.22 x 10~ 32 s 2 /s 3 f 4 l. Again, to 
be conservative, and to allow the filter to handle any possible instability resulting from clock 
“warm-up”, we set the initialization </2 value to 3.33 x 10 32 s 2 /s 3 . We kept the phase (bias) q 
unchanged. Below is a comparison of the old set and the two new sets of process noise values 
for Rubidiums: 

q Value OLD INITIAL NEW NORMAL 

Bias (</i) 1.11 x 10~ 22 s 2 /s 1.11 x 10" 22 s 2 /s 1.1 1 x 10~ 22 s 2 /s 

Drift (</ 2 ) 4.44 x 10“ 32 s 2 /s 3 3.33 x 10- 32 s 2 /s 3 2.22 x 10~ 32 s 2 /s 3 

Drift Rate (</ 3 ) 9.00 x 10 -42 s 2 /.s r ' 1.35 x 10~ 43 s 2 /s r ’ 6.66 x 10 _45 s 2 /s r ’ 

Of course, one might question using 30-day averaged drift rate values for deriving </ 3 — could 
the drift rate change by an unacceptable amount during those 30 days, thus undermining the 
premise of these calculations? Well, in the past, NRL has been able to apply as much as a 150 
day flat-average aging correction to their Allan Deviation plots-plots showing one-day stability 
figures similar to SVN25’sO°l. The implication is, if the Filter has a good drift rate term, that 
value can essentially be fixed for, in some cases, up to 150 days, without significantly degrading 
the one-day accuracy of the other clock states. Certainly, assuming drift rate consistency over 
30 days, let alone 150, was safe for deriving the above q values for the MC-S Kalman filter. 


SVN9 End Of Life Testing 

50th Space Wing approved a 1 SOPS and 2 SOPS joint effort to conduct End Of Life testing on 
SVN9 during March and April 1994 [19]. As part of the plan, Rockwell suggested dedicating 
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7-8 days for testing Rubidium clock drift rate estimation. We used the “New Normal” q values, 
and monitored the resulting system performance. 

The test, which lasted 8.7 days, produced very encouraging results! 19 !. At the end of the 
test, with tighter process noise values, the Kalman Filter converged on a drift rate value of 
—2.38 x 10~ 18 s/s 2 , with an associated standard deviation of 1.99 x 10 _19 s/s 2 (compared to a 
typical standard deviation of 1.0 x 10 -18 s/s 2 , using the old q values). Using an off-line tool, 
Rockwell derived post-processed values for comparison. Using a simple slope of their Aj 
(frequency) estimates over 7 days, Rockwell’s drift rate estimate was -2.44 x 10~ 18 s/s 2 , well 
within one sigma of the filter’s estimate. The National Institute of Standards and Technology 
(NIST) Report on SVN9 End of Life Testing pointed to a value of -2.32 x 10 -18 s/s 2 l*l, 
also well within one sigma of the Filter estimate. These comparisons indicated that the 
Filter had performed as designed-to converge on a more accurate drift rate estimate, with a 
correspondingly representative error estimate (standard deviation) [Figure 2], 

By, in effect, “clamping” on the filter estimate, one must question whether this covariance 
tightening is too restrictive, limiting the filter’s capability to respond to normal clock movement. 
We used two MCS parameters to test this capability. 


a. The first parameter was the Measurement Residual Statistical Consistency Test (MRSCT). 
Essentially, the MRSCT decides whether or not to accept Pseudoranges (PRs). Over 8.7 
days, the Filter accepted each and every smoothed PR for SVN9. The average PR residual 
(PRR) was no higher than that of a typical healthy, operational vehicle, or SVN9’s prior 
to the test. 

b. The second parameter was the Estimated Range Deviation (ERD). The ERD gives a 
good indication of the range error a user is experiencing, based on the current navigation 
upload residing in the vehicle. Over the 8.7 days, we uploaded SVN9 only once per day, 
and the ERD RMS never once exceeded 3.1 meters-well within our ERD criteria of 10 
meters. Correspondingly, the one-day User Range Accuracy (URA) dropped from 5.0 
to 3.8 meters, and the four-day URA dropped from 33.0 to 13.0 meters. In hindsight, 
we could have even set SVN9 healthy during the test, and netted a small improvement 
to global coverage and accuracy [Figure 2]. 


In short, results from the SVN9 drift rate test indicated that Filter estimation worked quite 
better with the reduced process noise (q) vaiues. 

Real World Implementation Of The New Rubidium qs 

On 18 Mar 94, we began applying these results towards real-world SVN10 and SVN24 clock 
estimation. Since, at that time, SVN24’s Rb w'as less than three months old, and since SVN10 
is a Block I, always susceptible to the effects of eclipse seasons, we selected the “Initialization” 
qs instead of the “New Normal” qs. 

For SVN10, during the three months prior to the test, ERDs exceeded 5.0 meters on 19 separate 
days. During the three months after the new qs were installed, SVN10 ERDs didn’t once 
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exceed 4.8 meters. In addition, our Smoothed Measurement Residual (SMRES) tool showed 
that SVN10 residuals from the DMA monitor stations, since 18 Mar 94, have been consistent 
with those prior to 18 Mar 94, as well as those for our other satellites. Similarly, between 
these two time periods, SVNIO’s time transfer error dropped from 14.6 to 9.9 nanoseconds 
(RMS), according to United States Naval Observatory (USNO) data! 5 '. These data points, 
from independent agencies, further show a significant improvement in satellite accuracy. 

Similar to SVNlO’s, the ERDs for SVN24 decreased after 18 Mar 94. Additionally, after 
installing the “New Normal” qs on 24 Apr 94, from that time to the present, the Filter has 
easily and consistently accepted SVN24 PRs. Likewise, SVN24 residuals from DMA, since 24 
Apr 94, have been as good or better than those prior to 24 Apr 94, and better than those 
of the other 23 operational satellites. In terms of upload accuracy, SVN24’s ERDs routinely 
exceeded 4.0 meters prior to 24 Apr 94. Since 24 Apr 94, SVN24’s ERDs have rarely exceeded 
3.5 meters, and have typically stayed under 2.5 meters. SVN24, now, is one of our two most 
accurate satellites. To complete the usefulness of this improvement, on 28 Apr 94, we included 
SVN24 into the GPS composite clock, allowing it to better stabilize GPS time. 

For the time being, after the 2 SOPS has initialized a Rubidium clock for 7-14 days, we’ll 
probably install the “Initialization” qs for 90 days. At the three month point, assuming nominal 
clock performance, we’ll likely install the “New Normal” qs. Also, at three months, we will 
aggressively consider including that satellite into the GPS composite clock — a Block II/IIA 
Rubidium clock estimate, now properly corrected for drift rate, now has a better one-day 
frequency stability than those of each of the on-orbit Cesiums. The GPS community, as a 
whole, can now at least tame a long existing ambivalence we’ve had about using Rubidium 
clocks in operational satellites. A Rubidium clock, now properly tuned in the Kalman Filter, 
significantly improves GPS timing and positioning accuracies. Currently, five GPS satellites use 
Rubidium clocks. One, in particular, SVN36, is arguably now our most accurate satellite. 

CESIUM CLOCK ESTIMATION 

Having resolved perhaps the most significant recent problem with GPS clock estimation through 
improved Rubidium qs, we decided to expand this opportunity for improvement to the remainder 
of all on-orbit GPS frequency standards: Cesium (Cs) clocks. As demonstrated earlier, deriving 
clock qs involves two main steps: 1) obtaining data that can accurately describe the behavior 
of the clocks involved, and 2) mathematically translating this behavior into the qs themselves. 

DMA has already been doing exactly this. A snapshot of some recently-derived DMA qs shows 
values that are, for the most part, unique to the individual clocks! 3 '. DMA’s qs vary significantly 
between satellites. In contrast, prior to 6 Oct 94, the MC.S qs were equal for most GPS Cs 
clocks. Also noteworthy is that the MCS’s qi value was less than each of DMA’s equivalent qj 
values! 3 !; 


67 



MCS q Values 


DMA q Values 


Smallest Largest 

Bias (</i) 1.11 x 10- 22 s 2 /s 4.25 x 10- 22 s 2 /s 4.35 x 10-' 21 s 2 /s 

Drift ( q 2 ) 4.44 x 10~ 32 s 2 /s 3 7.35 x 10 _33 s 2 /s 3 5.16 x 10 _32 s 2 /s 3 

This comparison raised two questions: 1) Would uniquely tuning the qs provide a significant 
improvement to GPS performance? 2) Does a legitimate reason exist for deliberately having 
lower qi terms in the MCS Kalman Filter? The remainder of this paper answers the first 
question. The second question, however, is more philosophical. 

MCS software experts will argue that a fundamental difference in purpose between the respective 
Kalman Filters at the MCS and at DMA constitutes a legitimate reason for using different 
q, terms. Since the MCS Kalman Filter is designed, in part, to provide accurate 24 hour 
predictions for navigation uploads, one could argue that we might want to deliberately keep 
our qi low to reduce the gain, and hence, prevent a situation whereby a noisy Kalman update 
could skew a 24-hour navigation upload prediction. Timing experts, however, will argue that 
tinkering with this parameter can be dangerous, since doing so can impose a configuration 
inconsistent with the basic intended design of a Kalman Filter. Both sides have very legitimate 
arguments. 

Deriving New Cesium Clock qs 

Analysts at NRL provide timely, accurate, and understandable reports on GPS clock perfor- 
mance. In particular, we now greatly utilize their Allan Deviation [ct 2 (t)] plots, created from 
DMA precise ephemeris data. The following equation relates the Allan Variance [<t 2 (t)] to 
Kalman Filter qs. This equation assumes independence between each sample frequency pair! 2 !. 

<t 2 (t) = </i (t _1 ) + <] 2 ( t )/3 + q3(?- 3 )/'20 (3) 


In order to relate current clock performance (via the Allan Deviation) to the system qs, we 
try not to use data more than 90 days old. Unfortunately, by only using 90 days of data, we 
experience the tradeoff of degraded confidence intervals for r > 20 days. For Cesium clocks, 
this is a non-concern, since we currently fix the drift rate and q 3 values to zero. For Rubidium 
clocks, however, the degraded confidence intervals, combined with the difficulty of correcting 
for drift rate without violating the sample frequency pair independence assumption, makes 
calculating the last term dangerous. As demonstrated earlier in this paper, we now have very 
suitable q 3 values for Rubidium clocks. Thus, for r < 20 days, we can substitute these into the 
Allan Variance equation, and simply solve for qi and q 2 - Then, we can compare our theoretical 
values to empirical values, using NRL Allan Deviation plots (with flat aging corrections applied 
for Rubidium clocks)( 12 l. 

One other concern relates to measurement noise. The data from NRL, and hence from DMA, 
has a fairly certain amount of measurement noise. The MCS’s parameter for measurement 


68 


noise, which we’ll call q 0 , accounts for some of the GPS monitor station (MS) receiver noise, 
some of the satellite clock’s white and flicker phase noise, MS location errors, and general 
modeling errors. DMA has a similar parameter designed to account for measurement noise, 
currently set to (45 cm) 2 S* 0.20 m 2 . For years, the MCS set this parameter at 1.0 m 2 . Thanks 
to recently refined MS location coordinates from DMA00, the MCS was recently able to reduce 
q 0 to (0.86 m) 2 = 0.74 m 2 . We derived this value using 500 Pseudorange Residual values from 
a widely distributed assortment of times and satellite-MS combinations. Our new value of 
0.74 m 2 , not surprisingly, is not a dramatic reduction from 1.0 in 2 , but nonetheless is consistent 
with our expectation of improvement from the new coordinates: 


y/(1.00 2 — 0.86 2 ) = .51 (meters) (4) 


One might suggest using DMA’s lower value. However, since our parameter accounts for more 
than just pure white measurement “noise”, our parameter is higher for a legitimate reason. 
Although not purely white phase noise in nature, noise associated with measurements can tend 
to misrepresent the stability of the estimated clock states. We can roughly express the instability 
resulting from this representation error asl 1 * : 


°'r( T ) = %>( T 2 ) 


( 5 ) 


By assuming independence between this representation error and the other noise processes on 
a given clock, the equation for the Allan Variance of the measured clock adds an additional 
term! 7 !: 


l7 m( r )=%( T 2 ) + </l( r ')/ :J + </3 (t 3 )/'20 


( 6 ) 


We created a Basic program to plot the theoretical a(r ) values, using the above equation, for 
r = 0.1 to 100 days. Using recent precise ephemeris rr(r) plots from NRU l2 l, along with the 
Basic program, we derived new q values for all satellites [Figure 3]. Note that the Rubidium 
qs remained unchanged. The Rb qs we derived earlier this year are, and have been, consistent 
with true clock performance. Nonetheless, Figure 4 shows how the theoretical Allan Deviation 
does change significantly for, in particular, SVN21 and SVN23, by using the newer qs. 

The current MS bias and drift qs, 1.11 x 10~ 22 s 2 /s and 4.44 x 10~ 32 s 2 /s 3 , respectively, are 
not representative of true MS clock performance. However, the MCS uses three separate 
mini-Kalman Filters, a.k.a. “partitions” to individually estimate MS clock states. Since a 
partition reconciliation algorithm keeps these states fairly consistent I Ml, over time, the MCS 
estimation structure effectively triples the weighting of the long term effects of MS clocks. With 
this current q 2 value for MSs, this "triple weighting” produces, in a roundabout fashion, the 
effect of using a q 2 roughly the same as the smallest satellite q 2 . We may tweak this parameter 
in the future, but, for the time being, this effect produces a fairly accurate resultf 16 '. 
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We also began using a newer set of qs during Cesium clock initialization. Below is a comparison 
of the old qs, and new initialization qs we’ve derived: 

Old ql Old q2 Old q3 New ql New q2 New q3 

10 -22 s2/s 10 _32 s 2 /s 3 10 -45 s 2 /s s 10- 22 s 2 /s l(r 32 s 2 /s 3 l(r 45 .s 2 /s r> 

1.11 4.44 0 4.44 3.33 0 


Testing The New Cesium q Values 

We safely tested the validity of these changes on 3 Oct 94, using a Test & Training simulator 
in the MCS. The results were impressive. 

a As expected, the state covariances converged to steady state values more truly represen- 
tative of the unique short- and long-term variances of the individual clocks. Also as 
expected, none of these new steady state covariances differed drastically from the typical 
older values. The implication of these small, but significant changes is that the Filter 
safely re-weighted clock state estimation based on true frequency standard performance, 
as opposed to assumed performance equality (equal qs): 

Value OLD VARIANCES NEW VARIANCES 

(All Cs) (Minimum) (Maximum) 

Bias 1.25 "x 10 _17 s 2 1.07 x 10“ 17 s 2 1.83 x 10 -I7 s 2 

Drift 3.20 x 10 _27 s 2 /s 2 1.38 x 10 _27 s 2 /s 2 4.38 x 10 _27 s 2 /s 2 

b As expected, the current state residuals experienced small (not trivial, not severe) changes, 
indicating that the Filter more responsibly distributed error to the appropriate states. 

c The MCS Pseudorange Residuals (PRRs) dropped from 1.61 m (RMS) to 0.87 m (RMS), 
after the Filter reprocessed the same raw data with the new set of qs. This more 
dramatically indicates that the Filter more responsibly distributed error to the appropriate 
states, so well that Filter predictions can now have less systematic error, and hence, less 
error when compared to smoothed measurements. 

d The consistency of MS clock states across the Kalman Filter partitions experienced a 
small, but not trivial improvement (A 3.8% reduction in Bias divergence error, and 21.6% 
reduction in Drift divergence error). Again, by more responsibly appropriating error 
to the respective clock states, short-term MS clock state instability across the partitions 
dropped. 


Real-World Implementation Of The New Cesium qs 

By installing these new qs on 6 Oct 94, we safely improved a) Kalman Filter clock estimation, 
b) navigation error representation, and c) the stability of the GPS composite clock. 
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The stability of GPS time, defined by the GPS composite clock, intuitively, should have improved 
simply as a result of the improved weighting, again, by uniquely tuning the qs based on true 
clock performance. When we used equal qs, the Allan Variance, <r^(r), of the implicit ensemble 
of N equally weighted clocks (for r = 1 day) was approximately! 4 !; 

4(r)-l/(A 2 )f;4(r) (7) 

i.= 1 

Using the one-day Allan Deviation figures from NRL Quarterly Report 94 — 3 r 15 1, the one-day 
stability of this implicit ensemble was approximately 1.55 x 10 -14 . 

By using clock-unique q s, the Allan Deviation of the now finely tuned implicit ensemble (for 
r = 1 day) is approximately! 1 ’ 4 !: 


i( T ) — 


£(*5i( T ))' 


l 


( 8 ) 


Incorporating the same one-day NRL Allan Deviation figures into the above equation, the 
one-day stability of the implicit ensemble dropped to approximately 1.22 x 10 -14 . Similarly, 
the observed Allan Deviation of GPS time, derived from USNO-smoothed measurements! 5 ’ 9 !, 
also dropped, not only for r = 1 day, but for 1 < r < 10 days [Figure 5]. 

Important to note is a large improvement in extended (14 day) navigation performance. By 
utilizing more representative (lower) q 2 values, the 14-day URA predictions have dropped to 
lower, more representative values for most satellites. Figure 6 shows a comparison of the 
typical 14-day URA values before and after 6 Oct 94, for all Block IIA satellites in estimating 
partitions. Though not an absolute indication of extended navigation accuracy, by uniquely 
tuning the qs, these URA values now, at least, have more validity than before. The 14-day URA 
values for all healthy GPS satellites, since 6 Oct 94, have been well below the NAVSTAR GPS 
System Operational Requirements Document (SORD) User Range Error (URE) specification 
of 200 meters! 17 !. 


CONCLUSION 

This fine tuning reinforces how deriving and installing clock-unique MCS Kalman Filter process 
noise values can safely and significantly improve GPS timing performance. We will continue 
to update these parameters on a regular basis. In the near future, we plan to review these 
values every three months, and as needed (after a clock swap or a dramatic change in clock 
performance). 

Loral Federal Systems Division received a tasking to more comprehensively review these and 
other data base parameters in 1995. We expect the results from their analysis to be more 
precise than the above results, due to the extensive background of the team of experts that will 
tackle this project. 
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Nonetheless, this successful attempt at fine tuning the MCS qs helps pave a path for future 
MCS data base analyses, and hence for future refinements to GPS timing performance. 
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QUESTIONS AND ANSWERS 


GERNOT M. WINKLER (USNO): Do you apply, or do you know whether anyone has 
every applied, any modeling for the ground stations, the same as you model your 12-hour 
satellite frequency change, which you observe? Has there been any modeling of the ground 
station? 

STEVEN HUTSELL (USAF): Ground station 12-hour periodics? 

GERNOT M. WINKLER (USNO): Whatever. 

STEVEN HUTSELL (USAF): Not that I’m aware of. I’m sure that there are seasonal and 
other effects. We try to model some of our tropospheric values. Right now, we have default 
values for those. LORAL Federal Systems Division has been working on that. The biggest 
problem we’ve had is getting environmental sensors that we’re comfortable with. Granted, we 
haven’t tested it as much as we could, and that’s probably something we could look into more. 
But for the time being, we are using database values. 
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